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ABSTRACT 

Following our initial discovery of blue luminescence in the spectrum of the Red Rectangle (RR) and its identi- 
fication as fluorescence by small three- to four-ringed polycyclic aromatic hydrocarbon (PAH) molecules, we 
report on the spatial correlation between the blue luminescence and the 3.3 fim emission, commonly attributed 
to small, neutral PAH molecules, and on the newly-derived UV/optical attenuation curve for the central source 
of the RR, HD 44179. Both results provide strong additional evidence for the presence of small PAH molecules 
with masses of less than 250 a.m.u. in the RR, which supports the attribution of the blue luminescence to fluo- 
rescence by the same molecules. We contrast the excellent spatial correlation of the two former emissions with 
the distinctly different spatial distribution of the extended red emission (ERE) and of the dust-scattered light 
within the RR. The UV/optical attenuation curve of the central star is unlike any interstellar extinction curve 
and is interpreted as resulting from circumstellar opacity alone. Major contributions to this opacity are absorp- 
tions in broad bands in the mid-UV, contributing to the electronic excitation of the luminescing PAH molecules, 
and a sharp ionization discontinuity near 7.5 eV in the far-UV, which places a sharp upper limit on the masses 
of the PAH molecules that are responsible for this absorption. The strength of the far-UV absorption leads to an 
abundance of the PAH molecules of 10~ 5 relative to hydrogen in the RR. Such small PAHs are perhaps unique 
to the environment in the RR, where they are shielded from harsh radiation by the dense circmstellar disk. 
Subject headings: ISM: individual(Red Rectangle) — ISM: molecules — dust, extinction — radiation mecha- 
nisms: general — stars: individual(HD 44179) 



1. INTRODUCTION 

Re cently, we repor ted the discovery of blue luminescence 
(BL ^Viih et al.ll2004T) in the spectrum o f the bipolar, proto- 
plane tary Red Rectangle (RR) nebula (ICohen et alJ 119751 
120041) . The emission, peaking near 378 nm, is contained in a 
band of FWHM ~ 45 nm. The band-integrated intensity of the 
BL is comparable to that of the scattered light intensity and of 
the intensity of the extended red emiss ion (ERE), which is 
also exceptionally stro ng in this object dSchmidt et al]ll980t 
IWitt & Borosonll990l) . This suggests that the BL is produced 
by a carrier of substantial abundance in the gas/dust medium 
in the RR. 

Through a comparison of the observed spectrum with lab- 
oratory photoluminescence spectra of a collection of likely 
dust components and with gas-phase fluorescence spectra of 
polycy clic aromatic hydrocarbon (PAH) molecules. lViih etafl 
J2004I) identified neutral three- and four-ring PAH molecules 
such as pyrene, anthracene and phenanthrene as the most 
likely sources of the BL. This identification also appears con- 
sistent with several other important facts. Firstly, the RR 
is the brightest known sou rce of emission in the unidenti- 
fied in frared (UIR) features jRussell et al.ll978tlGeballe et all 
which is commo nly attributed to PAH molecules and 
ions dBakes et alJ l200l which suggests high relative abun- 
dances and optimal excit ation conditions for such molecules. 
Also, as a post-AGB star llvan Winckel2003h . the carbon-rich 
central star of the RR, HD 44179, is in an active mass-losing 
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and dust-producing stage of its life which suggests that PAH 
molecules are currently co ndensing in the carbon-rich outflow 
llKellerfll987t ICherchneff et~aTlll992l This includes, in par- 
ticular, smaller PAH molecules that are able to exist in the 
relatively benign radiation environment of the RR, but which 
would not be viable after being ejected into the much harsher 
radiation field of interstellar space, wh ere three- and four- 
ringed PAHs are not expected to survive (tJochims et al.l!994l 
1999). And finally, anthracene and pyrene have exception- 
ally high flourescence efficiencies, which would cause their 
spectra to dominate any fluorescence spectrum of a mixture 
of PAH molecules jBrechignac & HermindH994T) . If con- 
firmed by further observations, the detection of three- and 
four-ringed PAH molecules would represent the detection of 
the largest specifically identified molecules so far observed 
outside the solar system. 

In this paper we report new observations and analysis deal- 
ing with the spatial distribution of the BL in the RR nebula. 
We show that the spatial distribution of the BL is distinctly 
different from that of the scattered light distribution in the 
same region of the spectrum and also distinctly different from 
that of the ERE. Significantly, there is an exceptionally close 
correlation between the distribution of the BL and the distri- 
bution of the UIR-band emission at 3.3 /im, attributed to the 
C-H stretch transition in neutral PAH molecules, which pro- 
vides strong support for the identification of the BL as fluo- 
rescence from relatively small PAH molecules. In the second 
part of this paper we discuss the newly-derived attenuation 
curve for the central star, HD 44179. The attenuation is shown 
to be circumstellar rather than interstellar in nature and to 
be dominated by absorptions attributable to electronic bound- 
bound and bound-free transitions in small PAH molecules. 
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2. OBSERVATIONS 

Low-resolution, long slit spectra of the RR were obtained 
with the R-C (Cassegrain) spectrograph at the Cerro Tololo 
Inter- American Observatory (CTIO) 1.5 m telescope. These 
observations were made on 2003 March 26 and 29, using a 
2'.' 5 wide slit, 7. '7 long. For observations in the blue (340 
- 600 nm), grating #09 with 300 1 mm" 1 was used which is 
blazed at 4000 A, provides a 8.6 A resolution and a spectral 
coverage of 260 nm. Using a CUSO4 filter to select the grat- 
ing's first order, the setup covered a wavelength range of 340- 
600 nm, including the range from H^ to the Balmer disconti- 
nuity. Grating #11 with 158 1 mm" 1 , blazed at 8000 A, with 
a 16.4 A resolution was used in the first order with a GG495 
cut-on filter for the red observations with a wavelength cover- 
age from 500-1000 nm. The Loral IK CCD detector yielded 
a spatial scale of 1"3 pixel" 1 along the slit. A coronographic 
decker assembly was used to minimize scattered light from 
the star while probing as much of the inner nebula as possi- 
ble. All observations were made using the full extent of the 
7. '7 long slit to get simultaneous sky observations. Spectra 
were taken at two nebular locations, 2."5 and 5" south of the 
central star HD 44179 with the slit in E-W direction (PA = 
90°). The nebular exposures were bracketed by exposures of 
the central star. Individual exposures were limited to 5-10 
minutes on the nebula and 1 minute for the star, and 4-5 ex- 
posures were obtained for each orientation. Data reductions 
were carried out with IRAF 2.12 EXPORT, and all spectra 
were flux calibrated via observation of sta ndard stars. 

The 3.3 /im data are reproduced from iKerr et alJ ( 119991) . 
The spectroscopic observations were made at the United 
Kingdom Infrared Telescope (UKIRT) on 1995 December 1 1 . 
The spectra were taken with a slit 90" long, 1 "2 wide, placed 
5" south of the central star HD 44179 (PA = -85°). 

For the FUV spectral energy distribution (SED) of HD 
44179 archival observations from the International Ultravio- 
let Explorer (IUE) LWP 22416, SWP 38188, and LWR 04273 
were used. 

David Malin (David Malin Images) provided us with one 
of the few existing deep blue images of the RR for the com- 
parison of our spatially resolved spectra with morphological 
details of the nebular structure. The image, reproduced in 
Figure 1, was taken by David Malin with the AAT at the 
Anglo- Australian Observatory. The exposure time was 15 
min, the detector was a Kodak Ila-O plate (blue sensitive) 
with a GG385 cut-on filter. The effective bandpass extends 
from 390 nm to 480 nm. A red image of the RR taken with 
the HST has been reproduced to aid in the correlation of the 
observed features with the nebular structure. 

3. ANALYSIS AND RESULTS 

3.1. Blue Luminescence in the Red Rectangle 

Nebular spectra were extracted in two-pixel bins along the 
two slits with 7 apertures along the 2 "5 south slit and 8 aper- 
tures along the 5" south slit with a spatial scale of 2"6 per 
aperture. Thus, each spectrum represents a 2 "5 x 2. "6 region 
of the nebula. Figure [2 shows the blue image of the nebula 
overlaid with the BL apertures and the 3.3 /xm slit and Fig- 
ure |2]is a red image overl aid with the central ERE apertures . 
The line-depth technique dViih et alJ20"04llWitt & Viihl2004l) 
was used to detect and measure the BL at the positions of 
the hydrogen Balmer lines for each spectrum. Measurement 
of the Balmer discontinuity in the nebular spectra compared 
to that in the stellar spectrum is used to infer the BL at the 




FIG. 1. — Blue image of the Red Rectangle with overlaid extraction aper- 
tures. The size of the image is 30" X 30". The white dot is an overlaid pixel 
to indicate the position of the central source. The dashed lines indicate the 
slit used for the 3.3 fim observation. North is to the top, east to the left. 



Balmer discontinuity (see Appendix A for details). Measure- 
ment of the equivalent widths of the absorption lines could 
also be used for a similar analysis, however a small uncer- 
tainty in the determination of the continuuium levels over the 
steep spectral shape results in a relatively large uncertainity in 
the equivalent width. In contrast the errors in the line-depth 
measurements are essentially limited to the photon noise in 
the spectrum. The co-added nebular spectra had a signal-to- 
noise ratio per pixel of ~ 10 in the outermost regions to ~ 
250 closer to the central star in the blue and a range of ~ 20 
- 225 in the red region of the spectrum. The stellar spectra 
had a signal-to-noise ratio per pixel of ~ 500 in both spectral 
ranges. 

3.1.1. Spectral Variability in the BL 

The BL spectrum is not identical in different regions of the 
nebula and certain trends can be noted. The pri mary peak in 
the spectrum is near 378 nm JViih et al.ll200l Fig. 4). A 
secondary peak around 397 nm starts to develop as we probe 
regions farther from the central star, more prominently along 
the 5" south slit. Also, farther from the star, the BL intensi- 
ties at longer wavelengths start to become stronger compared 
to the peak intensity. Interestingly, spectra that probe regions 
inside the conical outflow indicate that the BL in these regions 
has another peak at much shorter wavelengths, shortward of 
360 nm. Figure|5]depicts the BL spectrum at three such repre- 
sentative locations: 2"5 south, 7. "8 east, a position in the outer 
regions of the nebula; 2. "5 south, 2 "6 east, a region close to 
the central source, inside the bipolar outflow; 5" south, 5 "2 
west, a region along the farther slit just outside the cone wall. 

Based on the existing correlation between the wavelength 
of peak intensity of the fluo rescence spectru m and the molec- 
ular size of PAH molecules JViih etal.l2004h . these trends can 
be interpreted to suggest a size variation in the BL carrier at 
different locations in the nebula: closer to the central source 
we see an emergence of a smaller population of PAHs and far- 
ther out we see increasing contributions from larger emitters. 
Given the limited spectral resolution inherent in our detection 
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FIG. 2. — Red image of the RR with overlaid extration apertures. The 
size of the image is 19'.' 4 X 20"2. The HST image (Credits: NASA, ESA, 
Hans Van Winckel (Catholic University of Leuven, Belgium), and Martin 
Cohen (University of California, Berkeley)) has been reproduced to aid the 
identification of the nebular structure with the extracted apertures. North is 
to the top, east to the left. 
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FIG. 3. — Spectral variation of the BL. The solid line traces the intensity 
at 2'.'5 south, 7"8 east, the dashed line at 2"5 south, 2"6 east and the dash- 
dotted line at 5" south, 5"2 west positions. Solid circles indicate wavelength 
positions associated with the 2"5 south slit and solid squares indicate wave- 
length positions associated with the 5" south slit. Intensities are normalized 
to unity at the peak positions. 



technique and only fair spatial resolution these trends are at 
best an indication of the size distribution of the BL emitters. 
We defer a more complete analysis of the spectral variation of 
the BL spectrum until after completion of an ongoing program 
in collaboration with D.G. York which maps the complete RR 
nebula with dense spatial coverage using the 3.5 m telescope 
at Apache Point Observatory. 

3.1.2. Gradients of BL/Scatte red Light 

Figure shows the distribution of the ratio of band- 
integrated BL to the total scattered light in the same band at 
various positions along the two slits. Figure^ illustrates the 
spatial distribution of the BL and the scattered light intensi- 
ties on a normalized scale. Due to the anisotropy of the phase 
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function, scattered light is strongly forward directed and falls 
off steeply as we probe regions where our line of sight pene- 
trates regions of the nebula where larger scattering angles pre- 
dominate. The BL, on the other hand is due to fluorescence 
and is an isotropic emission and thus falls off less steeply. 
Thus, the spatial variation of the BL and scattered light along 
the two slits reveals the isotropic nature of the BL. There- 
fore, the ratio IsL/Iscat increases away from the star, whereas 
the actual intensities fall away from the star as the density of 
emitters and the exciting radiation both decrease. It is equally 
illuminating to note that the BL intensities have similar pro- 
files along the two slits, whereas the distribution of the scat- 
tered light it quite different; along the slit closer to the star, 
the scattered light falls off more steeply than along the second 
slit as the relative change in the distance from the central star 
varies more slowly in the latter case. 

3.1.3. Correlations with Other Emissions 

Figure [5] displays the spatial distribution of the band- 
integrated intensities of ERE and the BL along the two slits 
and the 3.3 fjm UIR-band emission along a 5" south (PA=- 
85°) slit (Fig. QJandElsh ow the apertu res). The 3.3 /im inten- 
sites (obtained from lKerr et alJ Jl999ft ~) have been normalized 
to the intensity of the BL at the east-west, 5"south posi- 
tion. This figure reveals a number of interesting facts. The 
3.3 /im UIR-band emission and the BL are almost perfectly 
correlated, which we will discuss in more detail in the next 
paragraph. The spatial distribution of the BL is distinctly dif- 
ferent from that of the ERE, the ERE is strongly peaked close 
to the central source and on the X-shaped walls of the outflow 
cone while the BL has a broader, more diffuse distribution. 

The 3.3 /im emission correlates exceptionally well with the 
BL along the 5" south slit, and the small differences are prob- 
ably attributable to the slightly different widths and PA of the 
two slits (see Fig. The 3.3 /im emission feature belongs 
to the distinctive set of mid-infrared emission features which 
have been attributed to emission from PAH molecules, which 
become vibrationally excited upon absorption of a UV pho- 
ton and subsequen t ly relax with the em ission of IR photons 
dBakes et alJl2004t ISchutte et aflll993l) . in addition to elec- 
tronic relaxation through fluorescence or phosphorescence 
transitions. Most of the vibrational transitions are insensi- 
tive to molecular sizes, but the 3.3 /im emission attributed 
to C-H stretch from PAH molecules is sensitive to the size 
of the molecular species. PAHs having ~ 20 carbon atoms 
show 3.3 lira emission intensities ~ 50 times stronger than 
that f rom PAHs having ^100 carbon atoms JSchutte et alJ 
Earlier observations of the RR showed a pronounced 
lack of similarity between th e 3.3 /im and the 11 .3 /im sur- 
face brightness distributions jBregman et aljfl993l) . the latter 
being attributed to C-H out-of-plane bending modes in rela- 
tively large PAH molecules. Fur thermore, the 3.3 emis- 
sion traces neutral PAH species jBakes et alJl200l , and be- 
comes extremely weak in ionized PAHs. The fact that the BL 
follows the 3.3 /im distribution so closely indicates that the 
BL emitters and the 3.3 /im emitters are similar small, neu- 
tral PAHs and thus strengthens our identification of the BL as 
fluorescence by small neutral PAHs. 

The ERE, BL and the 3.3 /im emission are all isotropic 
emissions and yet the spatial distributions of the ERE is dis- 
similar to those of the BL and the 3.3 /im emission. The ERE 
falls off much more steeply than the other two that have a 
broader distribution. The ERE is largely con fined to the bipo- 
lar, X-shaped structure of the nebula (see also lSchmidt & Witl 



Vijh et al. 



10 



0.1 



0.01 



2.5s-BL/Scat 
5s-BL/Scat 




10 



-15 



-10 



-5 

offset (arcsec) 



10 



15 



0.1 



0.01 



0.001 







2.5S-BL — •— 






5s-BL —■— . 






2.5s-Scat 






v\ 5s-Sc3t • — 










1 

















-15 



-10 



-5 

offset (arcsec) 



10 



15 



FIG. 4. — a. Ratio of band-integrated BL to scattered light along the two slits. Solid line traces the ratio along the 2'.'5 south slit and the dashed line that 
along the 5" south slit. b. Normalized BL and scattered light intensities. The BL is traced by solid lines and the scattered light by dashed lines. Intensities are 
normalized to unity at 2"6 east position. In both figures, solid circles denote spatial positions along the 2"5 south slit and the solid squares those along the 5" 
south slit. Positive and negative offsets indicate positions to the west and east respectively, along the slit. Zero offset refers to positions directly south of HD 
44179. 



H99ll Fig. 2,), whereas the BL is the dominant emission in 
the outer regions. The variation in the spatial distributions of 
these different emissions must be attributed the distribution 
of the carriers and the respective exciting radiations. The BL 
carriers, most likely small PAHs, are probably ionized inside 
the bipolar oufiow cone, where they are exposed to direct il- 
lumination from the central RR source. Thus fluorescence 
in these regions is effectively quenched. Farther out, outside 
the cones, the molecules/carriers are shielded from the far- 
UV ionizing radiation but still receive through scattering the 
mid-UV exciting radiation and thus emit with fairly high in- 
tensities. The ERE carriers on the other hand seem more ro- 
bust as the ERE intensities are high closer to the star inside 
the bipolar outflow cone and along the walls (X-shaped struc- 
ture). This could, in fact be interpreted as pointing toward 
some unknown ionized species as the carrier of the ERE. Out- 
side the outflow region the ERE intensities drop by over three 
orders of magnitude, suggesting that either the carriers of the 
ERE do not exist in such regions or if they do exist they do 
not receive the exciting radiation needed to produce the ERE. 

3.2. ATTENUATION OF HD 441 79 
3.2.1. Visual Reddening and Attenuation 

The central source of the RR, HD 44179, is a single- 
line spectroscopic binary star JWae lkens et aTl 119921 119961: 
van Winckel et al 1191)51 IMen'shchikov et alJ2002l) that is to- 
tally hidden from direct view at optical wavelengths by a 
circum-b inary disk seen nea rly edge on. Diffraction-limited 
imaging (ICohen et al.ll2004l) reveals the central source as a 
pair of scattered-light lobes above and below the disk, di- 
recting light from the hidden central stars via scattering by 
dust toward the Earth. The photometrically dominant com- 
ponent of the bina ry was classified as spectral type B9/A0 
llCohen et alJ 119751) originally, but is now recognized as a 
post-AGB star with an effective temperature near 8000 K 
with a high ly metal-deficient atmosphere typical of p ost-AGB 
stars (Ivan Winckel2003l iMen'shchikov et all2002l) . The low 
metal abundance, e.g. [Fe/H] = -3, gives this star the spec- 
troscopic appearance of a much hotter B9/A0 star, but with an 
intrinsic (B-V) color characterized by its actual effective tem- 



perature. By convolving the appropriate Kurucz atmosphere 
model (T eff = 8250 K, log g = 1.5, [Fe/H] = -3.0) with the 
B and V pass bands, we determine an intrinsic (B-V) = - 
0.04 mag for the photometrically dominant component of HD 
44179. We note that the nature of the secondary component 
is not known directly from observation, but its inferred mass 
of 0.35 M© and the fact that the innermost region of the RR 
contains a small H II region suggest that it is most likely a hot 
white d warf with a lumininosity of less than 2 % of that of the 
primary (|M en'shchikov et al 12001 . 

The observed (B-V) color of HD 441 79 has been inferred 
from broadband photometry of the RR bv lCohen et alJ (119751) 
as 0.39 mag., leading to an estimated color excess of E(B-V) 
= 0.35 mag. The peculiar nature of HD 44179, its loca- 
tion within an optically thick disk, and its uncertain distance 
make it difficult to estimate the total visual attenuation of its 
light as seen from Earth. Here we adopt the latest model 
(M en'shchikov et al 12001 which assigns a distance of 7 1 pc 
and a total luminosity of 6050 L Q to the central object, with 
uncertainties of 10 % a nd 20 %, respective ly. Given a visual 
brightness of V = 8.83 ( ICohen et alJl 19751) . we then estimate 
a visual attenuation of Ay = 4.2 mag. Here, we deliberately 
distinguish between attenuation and extinction. By extinction 
we understand the partial obscuration of direct star light due 
to absorption and scattering by dust along the direct line-of- 
sight to the star, while attenuation measures the reduction of 
observable flux due to absorption and scattering in a complex 
geometry, where the residual observable flux may indeed only 
be scattered light re-directed into the line of sight while the di- 
rect light is totally absorbed, as is the case in the RR. 

The attenuation of the light of HD 44179 is characterized 
by a ratio Ry = Ay/E(B-V) = 1 1.9, much larger than the av- 
erage value of Ry ~ 3.1 encountered with interstellar extinc- 
tion. The attenuation in the UV/visible portion of the spec- 
trum of HD 44179 is therefore much greyer than typical inter- 
stellar extinction, an effect attribut ed in the past to extin ction 
by larger grains in the RR ( IMen'shchikov et alJl200l . We 
consider it more likely that the large value of Ry arises from 
the fact that only scattered light from HD 44179 can be ob- 
served at optical wavelengths. Upon first scattering, scattered 
light is initially bluer than that of the illumination source, but 
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FIG. 5. — Spatial distribution of the BL and correlations with other emissions. Band-integrated intensities of BL and ERE at positions along the two slits (a. 
along the 2"5 south slit and b. along the 5" south slit) are plotted. ERE is plotted using dashed lines, the BL using solid lines and the 3.3 fim emission normalized 
to the value of the BL at 5" south, zero offset position, is plotted with a dotted line. Filled circles indicate positions along the 2"5 south slit and filled squares 
those along the 5" south slit. The intensities are band-integrated and are given in units of ergs cm~ 2 s~'sr~' for the BL and the ERE. 



upon transfer through an optically thick medium and multiple 
scattering, scattered light also reddens but at a slower rate than 
expected from the optical path length covered, as shown in the 
example of the colors of reflection nebulae by Witt (1985). 
The high value of the estimated visual attenuation of Ay = 
4.2 mag and the fact that scattered light must be re-directed 
by 90 degrees require that multiple scattering is the dominant 
process for the transfer of optical radiation from the center of 
the RR. 

3.2.2. UV/Optical Attenuation Curve 

In Figure[6]we display the observed (attenuated) UV/optical 
SED of HD 44179 over the wavelength interval from 130 nm 
to 600 nm. We produced this SED by combining the observed 
flux distributions in the near- and far-UV obtained from the 
IUE archive (LWP 22416, SWP 38188, and LWR 04273) with 
our optical SED obtained from CTIO spectroscopy. No ad- 
justments or corrections were applied; the continuity between 
the ground-based and space-based observations is almost per- 
fect, indicating that the two sets of absolute calibrations are 
in excellent accord. Also shown in Figure |6]is a mod el SED 
for a metal-deficient, giant-star using a iKuruczj (119931) stellar 
atmosphere for [Fe/H] = -3 , log g = 1 .5 and T e ff = 8250 K, 
which we propose as a suitable representation of the intrin- 
sic SED for the photometrically dominant component of HD 
44179. The model SED is normalized to the observed SED 
at V. The effective temperatu re of this model atmospher e is 
higher than that proposed by Men'shc hikov et all J2002I) by 
500 K, but our choice was dictated by the requirement that 
the model and the observed SED should have the same rel- 
ative strength of the Balmer discontinuity such that the ratio 
of these two SEDs would be continuous across the Balmer 
discontinuity. Model SEDs with temperatures differing from 
8250 K by as little as 250 K exhibit noticeable discontinuities, 
when divided into the observed SED. 

Figure displays the wavelength dependence of the 
UV/optical attenuation suffered by the light of HD 44179. 
The attenuation curve is normalized to E(B-V) = 1, and it 
is shown in comparison with the average Galactic interstel- 
lar extinction curve with identical normalization. Several dis- 
tinct differences between the two curves are immediately ap- 
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FIG. 6. — Observed UV/optical SED of HD 44179 compared to model 
SED. The observed SED has been obtained by combining the archived flux 
distributions in the near- and far-UV from IUE and optical flux distribution 
obtained by our CTIO spectroscopic measurements. The model SED is from 
a Kurucz (1993) stellar atmosphere for [Fe/H]=-3, log g =1.5, and T c ff=8250 
K, normalized to the obsereved SED at the V-band. The stellar fluxes are 
displayed in units of ergs citt's^'A -1 . 
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FIG. 7. — UV/optical attenuation curve for HD 44179. The average in- 
terstellar extinction curve is shown for comparison with the solid line. Also 
shown is the expected contribution of PAH ionization to the far-UV part of the 
attenuation curve by pyrene as a representative example with a dashed line. 
Indicated schematically are the four characteristic PAH absorption bands, the 
a, para, (3 and the f) 1 . The curves are normalized to E(B—V)=l. 



parent. With the slopes forced to coincide in the B-Vrange, 
the observed attenuation curve rises more steeply in the near- 
UV compared to the interstellar extinction curve. There is 
no sign of the 217.5 nm absorption peak; instead we see a 
very much broader attenuation hump reaching its peak near 
A -1 ~ 5 fim~ l . Finally, near A" 1 w 6.0 /im" 1 the attenuation 
curve exhibits a sharp discontinuous rise, which is much more 
abrupt than the gentle rise shown by the extinction curve. 

As is well known, the perception of anomalous attenua- 
tion/extincti on curves is stron gly dependent upon their nor- 
malization (tFitzpatrickl 12004ft . In Figure [8] therefore, we 
display the attenuation curve for HD 44179 normalized at 
V, scaled to Ay = 1.0, and compare it to Galactic extinc- 
tion curves with Ry = 3.1, representative of the diffuse ISM, 
and Ry = 5.5, representative of a molecular cloud envi- 
ronment and a grain size distribution dominated by larger 
grains. While the overall level of UV attenuation in HD 
44179 relative to the attenuation at longer wavelengths is 
roughly comparable to the Ry = 5.5 extinction curve and 
consistent with earlier findings of large dust grains in the 
RR dMen'shchikov et alJll998t lTura et alJll997l) . the two fea- 
tures that make the HD 44179 attenuation truly unique remain 
clearly visible: the absence of the 217.5 nm absorption peak 
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and the sharp discontinuous rise in attenuation at A" 
/im -1 . 

In two subsequent sections we will interpret the features 
just described in more detail. Given the overall appearance of 
the attenuation curve for HD 44179, we can conclude, how - 
ever, that the suggestion by iMen'shchikov et all J2002I) that 
most, if not all, of the reddening of HD 44179 is caused by 
interstellar instead of nebular material intrinsic to the RR is al- 
most certainly incorrect. The hydrogen column density (N(H 
I) < 1 x 10 cm" 2 ) implied by the equivalent width of the 
578 nm diffuse inte rstellar band in the spectrum of HD 44179 
JHobbs et aT]l2004l) translates to an extinction of A v < 0.05 
mag. Also noteworthy is the absence of narrow, interstellar 
absorption lines Na I and Ca II in the spectrum of HD 44179 
JHobbs et aT1l2004l) . The probability that the line-of-sight to 
one of the most unusual objects in the Galaxy should also 
be characterized by one of the most unusual interstellar ex- 
tinction curves ever observed is vanishingly small. Given the 
large internal optical dept h of th e RR implied by the visual 
attenuation found in Sect. 13.2.11 we suggest that the unique 
characteristics of the attenuation curve of HD 44179 must be 
explained in terms of the optical characteristics of the molec- 
ular and dust constituents of the RR. 

3.2.3. Far-UV Rise 
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The abrupt and discontinuous rise in the observed attenu- 
ation curve beginning at A" 1 pa 6 /Mm -1 , corresponding to a 
wavelength of 167 nm or a photon energy of 7.4 eV, is not 
characteristic of dust extinction but rather suggests the on- 
set of strong absorption, e.g. the absorption associated with 
the ionization discontinuit y in small PAH molecu les, which 
occurs at just this energy. IVerstraete et alJ (119901) published 
detailed ionization cross-sections for two PAH molecules, 
pyrene and coronene, which exhibit a sharp onset at A -1 pa 
6.0 /im" 1 , similar to the observed rise in the HD 44179 at- 
tenuation curve. When normalized to cross-section per C- 
atom, the absolute values of the ionization cross-sections and 
their wavelength dependencies are almost indistinguishable. 
A typical value is 1.35 x 10~ 18 cm 2 C-atom" 1 at 143 nm or 
A -1 « 7.0 ^m _1 . The work by iJochims et alJ (119961) shows 
similar ionization cross-section spectra for naphthalene, azu- 
lene, anthracene, phenanthrene, and benz(a)anthrathene, all 
with an sharp onset near 7.5 eV, corresponding to A -1 pa 
6.07 /LinT 1 . However, the ionization potential (IP) of PAH 
molecules is generally size-dependent in a manner similar 
to the size-dependence of the peak wavelength of fluores- 
cence (IViih et al.l 120041) . We demonstrate this dependence 
of the IP on molecular size by plotting known IP values 
( lEilfeld & Schmidtll981h against their molecular mass in Fig- 
ure|9l This graph illustrates that the rise in the far-UV attenua- 
tion curve is consistent with the ionization of PAHs, provided 
their molecular mass is not greater than about 250 a.m.u. A 
similar upper limit of the masses and siz es of PAH molecules 
was also deduced by IViih et aD J2004:) on the basis of the 
observed blue fluorescence seen in the RR. The attenuation 
curve suggests that absorption by molecules of this size and 
smaller could contribute to the sudden far-UV rise but that 
molecules of larger size cannot be present with significant 
abundance close to the central star; if they were, the disconti- 
nuity would occur at longer wavelengths. 

An important question is whether the observed increase in 
attenuation in the RR over the wavelength interval from A -1 pa 
6.0 /im -1 to A -1 ?a 7.0 /im -1 is in agreement with the likely 
carbon column density and the measured ionization cross- 
section at A" 1 pa 7.0 /zm" 1 . The observed increase in atten- 
uation over this interval is 2.0 mag , when reduced to the ac- 
tual color excess of E(B-V) = 0.35 of HD 44179, correspond- 
ing to an optical depth increase by At = 1.85. With an ion- 
ization cross-section a = 1.35 x 10~ 18 cm 2 C-atom -1 at 143 
nm wavelength, the corresponding carbon column density is 
Nc = 1.37 x 10 18 cm" 2 for carbon atoms tied up in small PAH 
molecules with masses less than 250 a.m.u. If the visual at- 
tenuation of 4.2 mag is associated with a mixture of gas and 
dust with relative elemental abundances as foun d in the Sun, 
then N (H I + 2H 2 )/A V = 1.87 x 10 21 cm" 2 mag- 1 ( iBohlin et al.l 
fT978b . and N(H I + 2H 2 ) = 7.85 x 10 21 cm" 2 . With a solar 
carbon abundance ratio of Nc/Nh = 3.4 x 10~ 4 , this leads to a 
total expected column density of carbon in all forms of Nc = 
2.7 x 10 18 cm" 2 . Hence, in order to produce the observed far- 
UV rise by ionization of small PAHs, we require that about 
50 % of all carbon atoms reside in small PAHs, if normal 
solar abundances prevail. Thus, the strength as well as the 
wavelength dependence of the observed far-UV rise in the 
attenuation curve are fully consistent with an origin in the 
ionization of small PAH molecules. In addition, HD 44179 
is recognized as a carbon-rich post-AGB star whose relative 
carbon ab undance must exceed that o f the Sun by factors of 
two to six (iCherchneff & Barkerll992l) . making the case even 
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FIG. 9. — Ionization potential of different families of PAHs as a func- 
tion of molecular mass. Also shown for r eference are the IP of graphite 
IClar & Schmidll973: lB urns & Yelke 1 969 ) and the ionization discontinuity 
apparent in the attenuation curve of the central source of the RR. 



stronger. In Figure we have indicated the expected contri- 
bution of PAH ionization to the far-UV part of the attenuation 
curve for pyrene as an representative example. The fact that 
the observed rise is slightly steeper than the individual curve 
of a specific PAH molecule can be understood easily, if still 
smaller PAH molecules with higher IP values contribute to the 
absorption at the shorter wavelengths. 

3.2.4. Mid-UVHump 

Given the facts that blue fluorescence from PAH molecules 
with an intensity comparable to that of the ERE contributes 
strongly to the the nebular emission of the RR and that the 
ionization of these molecules in the inner part of the RR 
contributes strongly to the far-UV attenuation at A -1 > 6 
fim~ l , it appears most likely to us that the mid-UV hump 
in the attenuation curve is produced by absorption related to 
electronic bound-bound transitions in these same molecules. 
With a sharp onset of the fluorescence near 360 nm, the ab- 
sorption powering this emission m ust occur at shorter wave- 
lengths. According to lClarl ( 119721) . the absorption spectra of 
PAH molecules consist of four characteristic bands, the a, 
para, f3, and j3' bands, located at progressively shorter wave- 
lengths in the mid-UV. They are schematically indicated in 
the mid-UV portion of Figure with arrows starting at the 
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band edges and the height indicating the relative strengths of 
the bands (See also ISchmidtll 1 977T) . The /3-band which typ- 
ically occurs at a wavelength 1/1.35 of that of the a band, 
is typically two orders of magnitude stronger than the a- 
band, while the intermediate para-band is typically one or- 
der of magnitude stronger than the a-band. The wavelength 
positions of these bands are closely related to the IPs of th e 
corresponding molecules dClar & Schmidl 1975lfl977l.il 9791) . 
so that a continuum of sizes of small PAH molecules with 
IP > 7.4 eV can be expected to produce a broad, quasi- 
continuous absorptio n band as observed in the attenuation 
curve of HD 44179. IVerstraete & Leg erl ill 9921) produced a 
generic mid- and far-UV absorption spectrum for interstellar 
PAH molecules, which exhibits the mid-UV hump and the 
far-UV rise so clearly apparent in the attenuation curve of HD 
44179. 

4. DISCUSSION 
4.1. Identification of the BL Carrier 

When we first attributed the BL in the RR to fl uorescence 
by thr ee- and four-ringed neutral PAH molecules dViih et alJ 
120041) . the identification of the likely carrier was tentative for 
several reasons. Our observational technique for separating 
the BL spectrum from the underlying scattered light spectrum 
produces measurements of the BL intensities only at the po- 
sitions of the hydrogen Balmer lines, yielding data with very 
limited and uneven spectral resolution that is incapable of re- 
vealing fine spectral features needed for specific identifica- 
tions of individual molecular species. In addition, the fluo- 
rescence spectra of several likely molecular PAH species oc- 
cupy nearly the same wavelength interval, contributing to the 
observed spectrum in proportion to their relative abundances 
and their fluorescence quantum yields. Finally, the gas-phase 
fluorescence spectra of PAHs vary as a function of temper- 
ature, with s ub-bands losing their sharpness with increasing 
temperature l lChi et al1l2001l) . For these reasons, the identifi- 
cation of the carrier of the BL can only be made in part on the 
basis of spectral comparison and must rely largely upon other 
supporting evidence such as presented in this paper. 

There are prospects for improving the spectral resolution 
of the BL spectrum by using other techniques for separat- 
ing the BL spectrum from the scattered light spectrum. One 
such technique is spectropolarimetry. It relies on the fact 
that the scattered n ebular light is highly linearly polarized 
JPerkins et al.lfl98l while the BL is expected to be unpo- 
larized. The presence of the BL in the composite spectrum 
of the RR would then be revealed by a proportional reduc- 
tion in the linear polarization of the combined light. The 
only two published sets of spectropola rimetric observations 
at optical wavelength s of the RR nebula JSchmidt et alJl980l 
iReese & Sitkdfl996l) do indeed reveal a sudden reduction of 
the linear polarization shortward of 410 nm consistent with 
the onset of the BL band. Unfortunately, the existing data 
lack the spectral resolution and adequate spectral coverage at 
shorter wavelengths that would be necessary. A program of 
spectropolarimetric observations of the RR is being planned. 

At this stage, the strongest supporting evidence for the cor- 
rectness of our identification of the carrier of the BL are the 
close spatial correlation with the 3.3 /xm C-H stretch emis- 
sion that is highly specific to small, neutral PAH molecules, 
and the evidence obtained from the analysis of the attenuation 
curve of the central source that exhibits strong features identi- 
fiable with electronic-band absorptions and an ionization dis- 
continuity specifically pointing toward PAH molecules with 



masses of 250 a.m.u. or less. 

4.2. Spectral Variability of the BL 

Our initial results showing evidence of systematic spectral 
variability of the BL with position within the bi-polar struc- 
ture of the RR are intriguing but clearly need more observa- 
tions with the highest possible spatial and spectral resolution. 
Such observations would have the potential of probing the 
variation in the dominant sizes of the PAH molecules as a 
function of distance from the central source and as a function 
of position within and outside of the bi-polar outflow cones 
and within the shadow of the circum-source disk. Our initial 
results are consistent with an increase in the size of the most 
strongly fluorescing PAH molecules with increasing distance 
from the central source, with a superimposed azimuthal de- 
pendence with respect to the outflow axis of the RR. This 
could reflect the growth of the molecules with increasing 
distance, but it could also be a result of the different exci- 
tation requirements for the fluorescence of PAH molecules 
of different sizes. The wavelengths of the absorption bands 
of PAH m olecules are related to their r espect ive ionization 
potentials dClar & Schmidtlll975L 119771 119791) and occur at 
shorter wavelengths for smaller PAHs, and at longer wave- 
lengths for larger PAHs. As a consequence, exciting radiation 
for the smallest PAHs may be lacking at increased distances 
as a result of the high opacity of the intervening nebular ma- 
terial. 

4.3. Spatial Variation of the BL 

The close spatial correlation between the BL and the 3.3 
fjm emission along two nearly identically placed slits, posi- 
tioned 5"south of HD 44179 is remarkable in itself by sup- 
porting the interpretation that both emissions are produced by 
the same molecules. It is further remarkable in that both emis- 
sions show a distinctly different distribution than the distribu- 
tion of the scattered light in the RR. The latter is produced by 
dust with a strongly forward-directed scattering phase func- 
tion, which when coupled with an embedded source, produces 
a strongly peaked brightness distribution of the resulting scat- 
tered light. This is clearly apparent when viewing the blue 
image of the RR (Fig. Q, which is dominated by a bright 
circular source of about 8" diameter containing most of the 
scattered light. The PAH fluorescence, both electronic in the 
form of the BL as well as vibrational in the form of the 3.3 
fjm emission, by contrast is much less centrally peaked, con- 
sistent with an isotropic emission process for both. It is fur- 
ther remarkable that the BL brightness distribution does not 
reflect the crossing of the X-shaped arms or wh iskers of the 
RR, w hich are so prominent at red wavelengths l lCohen et alJ 
120041) . This indicates that the BL is neither associated with the 
interior of the outflow cones nor with the walls of the outflow 
cavity but rather is associated with a population of emitters re- 
siding in a shell surrounding the inner bi-polar structure. This 
would place them into the shadow of the circum-source disk 
and thus outside the reach of ionizing radiation, consistent 
with our attribution of the BL to neutral PAHs. 

The difference between the spatial distributions of the BL 
and the ERE is striking. It has been known for some time that 
the ERE is strongly concen trated in the X-shaped arms of the 
RR ISchmidt & Wittlll99lh . which are thought to be the pro- 
jections of the walls of the bi-polar outflow cavity upon the 
plane of the sky. Our data leave little doubt that the two emit- 
ters involved require totally different excitation conditions. 
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Given the ubiquity of the ER E in a wide range of interstel- 
lar environments (IWitt & Viihll2004l) . it is not likely that the 
concentration of the ERE in the walls of the RR outflow cav- 
ity is the result of an abundance enhancement but rather re- 
flects the specific excitation re quirement of the ERE process. 
As reviewed by JWitt & Viihl 120041) . there are several inde- 
pendent indicators suggesting that photon with energies E > 
7.25 eV are required to initiate the ERE. These energies are 
much higher than those required for the excitation of fluo- 
rescence in neutral PAH molecules and the number of such 
high-energy photons are much more limited in the SED of 
HD 44179. Thus, they can barely penetrate the walls of the 
outflow cavity before their number is exhausted, thus limiting 
the spatial extent of the ERE-producing emitters. The topic of 
the ERE excitation will be the subject of another publication 
by the present authors. 

4.4. Attenuation of HD 44179 

The attenuation curve of HD 44179 appears highly unusual 
when compared to the wavelength dependence of interstel- 
lar extinction, but it appears to be fully consistent both qual- 
itatively and quantitatively with expectations, if the scattered 
light of HD 44179 has to penetrate a large column density (Ay 
= 4.2 mag) of circumstellar material filled with neutral PAH 
molecules with atomic masses less than 250 a.m.u. 

With a column density of carbon in the form of small PAH 
molecules estimated from the strength of the far-UV absorp- 
tion related to PAH ionization, we can estimate a molecular 
abundance of such PAH molecules with typically 16 carbon 
atoms per molecule relative to hydrogen in the RR of about 
10~ 5 . This agrees well with the canonical abundances of PAH 
molecules derived from the analysis of the strength of the mid- 
IR aromatic emission features, com monly referred to as the 
unide ntified infrared emission bands tLeger & d'Hendecourtl 
qui. This makes PAH molecules the most abundant species 
of interstellar molecules after H2 and CO. The fact that despite 
the strong indications for a presence of small PAH molecules 
there is no clear counterpart for the interstellar 217.5 nm ab- 
sorption feature in the attenuation curve of HD 44179, sug- 
gests that proposals to associate the 217.5 nm feature with 
PAH absorption in current dust models should be viewed with 
caution. The current observations appear to prohibit at least 
a role for small PAHs with masses < 250 a.m.u. This still 
leaves the possibility that larger PAH structures could be in- 
volved in contributing to the interstellar 217.5 nm feature, 
because they would have the advantage of greater thermody- 
namic stability under interstellar conditions. The small PAH 
molecules apparent in the RR are able to form and survive in 
the relatively benign radiation environment of the RR, espe- 
cially in the shadow of the optically thick dust ring blocking 
direct radiation from the central source. It is in these parts of 
the RR nebula where the BL intensity is strongly enhanced 
in comparison to the scattered radiation and the ERE. These 
small PAH molecules, however, are not likely to survive the 
harsher radi ation conditions prevailing in the diffuse interstel- 
lar medium JJochims et all 199411 19991) . where they would be 
subject to photo-dissociation. As a consequence, it does not 
appear likely that BL with spectral characteristics and intensi- 
ties similar to those observed in the RR will be seen in dusty 
star-fo rming regions such as reflection nebulae and H II re- 
gions jRush & Will 19751) . 

The conditions deduced for the RR environment also places 
interesting constraints on the carriers of the diffuse interstellar 
bands (DIBs). Current ideas for the carriers of these much- 



studied but still unidentified interstellar absorption features 
concentr ate on carbon aceous molecules, their ions in partic- 
ular (see ISnowlll995l for a review). PAH ions have been 
considered as the most likely carriers by a number of inves- 
tigators. Given the strong indicators for a large column den- 
sity of PAH molecules and their ions in the RR, it is remark- 
able that on ly one DIB A5780 is detect e d in th e spectrum of 
HD 44179 (iHobbs et al.ll200H) . ISnowl (12001 has proposed 
that this fact might be explained, if DIBs are produced by 
PAH cations with masses well in excess of 250 a.m.u. An- 
other fact that must be taken into consideration is the ab- 
sence of absorptions resulting from the D2 ^Dq transition of 
anthr acene* (708.76±0 13 nm ), naphthalene + (670.73±0.06 
nm) dSuk horukov et al.| 20041) . and acenaphthene" 1 " (646.3 
rim) llBieTmTer"etldjl27)03 ) which might be expected in the 
spectrum of HD 44179, if these cations are present in the light 
path. Some of these discrepancies could be resolved by con- 
sidering the ionization of these PAHs to the dication stage. 

5. CONCLUSIONS 

We summarize our conclusions: 

1. Our spatially resolved observations of the BL in the 
RR have yielded the first tentative evidence that the 
detailed spectrum of the BL is spatially variable, con- 
sistent with a change in the dominant size of the fluo- 
rescing molecules with position relative to the central 
source. 

2. The spatial distribution of the band-integrated BL in 
the RR differs fundamentally from those of the dust- 
scattered radiation and of the ERE, consistent with an 
origin from a nebular volume located mainly outside 
the outflow cavity and concentrated in the shadow of 
the disk obscuring the central source. This supports our 
earlier conclusion that the BL is produced by fluores- 
cence from small, neutral PAH molecules. 

3. There exists an excellent spatial correlation between the 
BL and the distribution of 3.3 /iin C-H stretch emission 
from PAHs. The latter emission is predominantly pro- 
duced by small, neutral PAH molecules, and the close 
spatial correlation supports the suggestion that the BL 
has its origin in the similar molecules. 

4. By combining our observations with existing IUE data 
for HD 44179, we have determined the UV/optical at- 
tenuation curve for the central source in the RR. The 
attenuation curve is characterized by an exceptionally 
high value of Ry = AyfE(B-V) = 1 1 .9, a remarkable ab- 
sence of the 217.5 nm absorption bump, the latter being 
a characteristic of Galactic interstellar extinction, and a 
strong, discontinuous rise in attenuation near A -1 « 6.0 
/im . These unusual characteristics suggest strongly 
that the reddening observed in the light from HD 44179 
is a result of radiative transfer in the RR and not due to 
interstellar extinction. 

5. The far-UV rise in the attenuation curve of HD 44179 
is qualitatively and quantitatively different from the 
much more gradual far-UV rise in Galactic extinction 
curves. The rise is, however, consistent with the on- 
set of photo-ionization of small PAH molecules with 
masses of less than 250 a.m.u., if the abundance of such 
PAH molecules relative to hydrogen is about 10~ 5 . This 
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represents a fully independent confirmation of the pres- 
ence of small PAH molecules in the RR environment, 
consistent with our conclusion regarding the origin of 
the BL in the RR. 

6. The broad mid-UV hump in the attenuation curve of 
HD 44179 whose appearance is quite unlike that of the 
familiar 217.5 nm interstellar absorption band, can be 
explained satisfactorily as resulting from a superposi- 
tion of the (typically) four mid-UV absorption bands 
characteristic of neutral PAH molecules. An observ- 
able consequence of this absorption and the resulting 
electronic excitation is de-excitation via fluorescence, 
which we indeed observe. 

7. We do not expect that BL with comparable relative in- 
tensities as seen in the RR will be seen in star-forming 
regions such as reflection nebulae and HII regions, 
where the mid-IR UIR band emissions are otherwise 
strong. The small PAH molecules required by the BL 
are not expected to survive in the radiation fields typical 
of star-forming regions. 
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APPENDIX 



DETERMINATION OF BL INTENSITY AT THE BALMER IUMP 



Assume that the ratio of the scattered light intensity at A less than the Balmer jump (BJ) to the scattered light intensity at A 
greater than the BJ, /<bV^>bj is determined by the stellar spectrum. We measure A<bj at 357 nm and A>bj at 397 nm. Assuming 



a constant albedo over the wavelength ran ge and an optically thin scenario, the relative increase in the scattered flux at 357 nm 
compared to that at 397 nm is 1.08 ( Savage & Mathis 1979). This factor is an upper limit as it will be smaller in an optically 
thick case. Thus, in the stellar spectrum 



TSC* 

'357 
'397 



: 0.2334 



After correction for wavelength dependent scattering, 



JSC* 

'357 
'397 



= 0.2520 



I$ = 0.2520*/^ 
where the superscript TV denotes the values for the nebular spectra. Also, 



Therefore, the intensity of the BL at 357 nm, 



jtotaIN 
'357 " 



jBL , jscN 
'357 +'357 



= + 0.2520 



jscN 



jBL jtotaIN tscN , r> oc-»n 
'357 = '357 _ '397 * 0.2520 



which can be determined by knowing the value of the scattered light intensity at 397 nm, which in turn is previously determined 
by subtracting the BL intensity (measured using the line-depth technique) at that wavelength from the total intensity at that 
wavelength. 



REFERENCES 



Bakes, E.L.O., Bauschlicher, C, & Tielens, A.G.G.M. 2004, in Astrophysics 

of Dust, ASP Conf. Ser. 309, ed. A.N. Witt, G.C. Clayton, & B.T. Draine, 

(San Francisco: ASP), 731 
Biennier, L., Salama, F, Allamndola, L. J., & Scherer, J. J. 2003, J. Chem. 

Phys., 118,7863 
Bohlin, R.C., Savage, B.D., & Drake, J.F 1978, ApJ, 224, 132 
Breginac, P., & Hermine, P. 1994, in AIP Conf. Proc. 312, Molecules and 

Grains in Space, ed. I. Nenner (New York: AIP), 613 
Bregman, J.D., Rank, D., Temi, P., Hudgins, D., & Kay, L. 1993, ApJ, 411, 

794 

Burns, J., & Yelke, E. 1969, Rev. Sci. Instrum., 40, 1236 

Cherchneff, I„ & Barker, J. B. 1992, ApJ, 394, 703 

Cherchneff, I., Barker, J.B., & Tielens, A. G. G. M. 1992, ApJ, 401, 269 



Chi, Z., Cullum, B. M., Stokes, D.L., Mobley, J., Miller, G.H., Hajaligol, 

M.R., & Vo-Dinh, T. 2001, Spectrochem. Acta A, 57, 1377 
Clar, E., & Schmidt, W. 1975, Tetrahedron, 31, 2263 
Clar, E., & Schmidt, W. 1977, Tetrahedron, 33, 2093 
Clar, E., & Schmidt, W. 1979, Tetrahedron, 35, 1027 

Clar, E. 1972, The Aromatic Sextet, London: John Wiley & Sons, Chapter 4 
Cohen et al. 1975, ApJ, 196, 179 

Cohen, M., van Winckel, H., Bond, H.E., & Gull, T.R. 2004, AJ, 127, 2362 
Eilfeld, P., & Schmidt, W. 1981, J. Electron Spectrosc. Rel. Phen., 351, 101 
Fitzpatrick, E.L. in Astrophysics of Dust, ed. A.N. Witt, G.C. Clayton, & B.T. 

Draine, ASP Conf. Ser. 309, 2004 (San Francisco: ASP), 33 
Geballe, T.R., Lacy, J.H., Person, E.E., McGregor, P.J., & Soifer, B.T. 1985, 

ApJ, 292, 500 



Small PAHs in the Red Rectangle 



11 



Hobbs, L. M. et al. 2004, ApJ, in press 

Jochims, H. W., Baumgaertel, H., & Leach, S. 1996, A&A, 314, 1003 
Jochims, H.W., Ruhl, E., Baumgaertel, H., Tobita, S., & Leach, S. 1994, ApJ, 
420, 307 

Jochims, H.W., Baumgaertel, H., & Leach, S. 1999, ApJ, 512, 500 
Jura, M., Turner, J., & Balm, S.P. 1997, ApJ, 474, 741 

Keller, R. 1987, in Polycyclic Aromatic Hydrocarbons and Astrophysics, 
NATO ASI Ser. C191, ed. A. Leger, L. d'Hendecourt, & N. Boccara, 
(Dordrecht: Reidel Publ. Comp.), 387 

Kerr, T.H., Hurst, M.E., Miles, J.R., & Sarre, P.J. 1999, MNRAS, 303, 446 

K urucz, R. 1993 CD-ROM # 13, CDBS from STScI. 
http : / /www . st sci . edu/hst / observatory /cclbs /k 9 3models 

Leger, A., & d'Hendecourt, L. 1987, in Polycyclic Aromatic Hydrocarbons 
and Astrophysics, eds. A. Leger, L. d'Hendecourt, & N. Boccara, 
(Dordrecht: Reidel Publ. Comp.), NATO ASI Ser., p. 223 

Men'shchikov, A B., Balega, Y.Y., Osterbart, R„ & Weigelt, G. 1998, New 
Astron., 3, 601 

Men'shchikov, A.B., Schertl, D., Tuthill, P.G., Weigelt, G., & Yungelson, L.R. 

2002, A&A, 393, 867 
Perkins, H.aG., Scarrott,5S.aM., Murdin,SP, & Bingham.aR.SG. 1981, 

MNRAS, 196, 635 
Reese, M.D., & Sitko, M.L. 1996, ApJ, 467, L105 
Rush, W.F., & Witt, A.N. 1975, AJ, 80, 31 
Russell, R.W., Soifer, B.T., & Willner, S.P. 1978, ApJ, 220, 568 
Savage, B.D., & Mathis, J.S. 1979, ARA&A, 17, 73 
Schmidt, W. 1977, J. Chem. Phys., 66, 828 



Schmidt, G.D., Cohen, M., & Margon, B. 1980, ApJ, 239, L133 
Schmidt, G.D., & Witt, A.N. 1991, ApJ, 383, 698 

Schutte, W. A., Tielens, A. G. G. M., & Allamandola, L. L. 1993, ApJ, 415, 
397 

Snow, TP. 1995, in The Diffuse Interstellar Bands, ed. A.G.G.M. Tielens and 
TP. Snow, (Dordrecht: Kluwer Academic Publishers), 379 

Snow, TP. 2004, Presentation at the 59th Symposium on Molecular 
Spectroscopy, Columbus OH, 24 June, 2004 

Sukhorukov, O. et al. 2004, Chem. Phys. Lett., 386, 259 

van Winckel, H. 2003, ARA&A, 41, 391 

van Winckel, H., Waelkens, C, & Waters, L.B.F.M. 1995, A&A, 293, L25 
htMej straete, L, & Leger, A. 1992, A&A, 266, 513 

Verstraete, L., Leger, A., d'Hendecourt, L., Dutuit, O., & Defourneau, D. 

1990, A&A, 237, 436 
Vijh, U.P., Witt, A.N., & Gordon, K.D. 2004, ApJ, 606, L68 
Waelkens, C, van Winckel, H., Trams, N.R., & Waters, L.B.F.M. 1992, A&A, 
256, L15 

Waelkens, C, van Winckel, H., Waters, L.B.F.M., & Bakker, E.J., 1996, 

A&A, 314, L17 
Witt, A.N. 1985, ApJ, 294, 216 
Witt, A.N., & Boroson.T.A. 1990, ApJ, 355, 182 

Witt, A.N., & Vijh, U.P. 2004, in ASP Conf. Ser. 309, Astrophysics of Dust, 
ed. A.N. Witt, G.C. Clayton, & B.T Draine (San Francisco: ASP), 115 



